FLOOD STUDIES IN TWO CHANNELS 

WITH A LINK 


A Thesis Submitted 

in Partial Fulfillment of the Requirements 
for the Degree of 

Master of Technology 



to the 

DEPARTMENT OF CIVIL ENGINEERING 
INDIAN INSTITUTE OF TECHNOLOGY KANPUR 


May, 2005 



*'8 'jui, 

uroftf 1 WR •wijpi 

Wl#i Vo A...I.lilEEIL^«<» 



CERTIFICATE 


It is certified that the work contained in the thesis entitled “Flood Studies in Two 
Channels with a Link”, by Kedar sharma, has been carried out under my supervision 
and this work has not been submitted elsewhere for the award of a degree. 



(Dr. Rajesh Srivastava) 
Associate Professor 


Department of Civil Engineering 
Indian Institute of Technology Kanpur 
Kanpur, India 



(Dr. P.K.Mohapatra) 

Assistant Professor 
Department of Civil Engineering 
Indian Institute of Technology Kanpur 
Kanpur, India 


Dated: May, 2005 




ABSTRACT 


Analysis of flood wave propagation is an important engineering practice. In this work, 
flood wave propagation is studied in two channels with a link. For the purpose, one- 
dimensional equations for water and sediment flow are considered. A high-order shock- 
capturing Total Variation Diminishing (TVD) scheme is used for modeling water flow. In 
addition, the sediment equation is solved by using a different scheme with a semi-coupled 
approach. This model can simulate either steady or unsteady flows. In case of rigid bed 
channels, the numerical model is validated against earlier analytical results (Stoker 
solution) for dam-break flow. However, the present model is validated against 
experimental results for sediment overloading, in case of mobile channels. The numerical 
model is applied to two channels with a link to simulate the unsteady flow in both 
channels. In the absence of experimental data, the results are compared with those 
obtained from a popular computer package (HEC-RAS). In all the cases, the match 
between the results is satisfactory and the present model captures the sharp wave front in 
a better manner. The numerical model is also used to simulate flow distributions in case 


of flows in rigid bed channels. 
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CHAPTER I 

INTRODUCTION 


All human civilizations have been associated with rivers; all major cities in the world are 
situated on rivers. Therefore, flood studies occupy an important place in engineering 
application. If the geographical conditions are favorable, then two channels can be linked 
to transfer the water from one channel to the other. Linking of two channels may help in 
flood control and in surplus water distribution. In an irrigation system, interlinking of 
canals is very common. 

Many schemes of large scale water transfer projects (interlink proposals) have been 
planned; some of them are implemented and are standing as a land-mark for the over- all 
development of the respective regions. In this transfer, water may flow by gravity or by 
pumping. In Canada, sixteen inter-basin water transfer schemes have been implemented 
for hydropower development. In USA, the longest and best known schemes implemented 
so far is California State Water Project, which envisages transfer of water from 
Sacramento river in North California to southwards through a 715 km long aqueduct with 
a lift of about 1000 m to meet domestic, industrial and irrigation demands. Similarly, in 
China, there are schemes existing from ancient times and recently, supplemented by 
modem construction techniques. China is also planning for transfer of 48 BCM of water 
from South to North through Grand Canal, close to the Eastern Coast. In India, many 
interlinking projects were implemented for hydropower and irrigation development. 
Periyar Project, Parambikulam Aliyar, Kumool Cudappah Canal, Telugu Ganga Project, 
Ravi-Beas-Sutlej- Indira Gandhi Nahar Project are some of them. The Sardar Sarovar 
Project on the river Narmada in Gujarat has successfully linked ten rivers of the state and 
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has provided a roadmap for the national river-interlinking project. The Narmada Main 
Canal (NMC), with a capacity to carry 40,000 cusecs of water, has played a major role in 
channelising the surplus water from the river Narmada to other deficit river basins. A brief 
note about the interlinking of rivers in India is given in appendix. 

If the two channels are linked, the flow is affected in both the channels. The water 
transferred from one channel may create flood condition in the other. It is also possible 
that a link makes the water deficiency in the donor channel. After the linking of channels, 
design of any engineering structure on these two channels is done by taking into account 
the effect of link on flow. The link also affects the sediment and pollutant transport in the 
channels. 

This problem can be studied by experimental, analytical or numerical modeling (Tannehill 
et al. 1997). Experimental approach is more realistic than tuialytical and numerical 
modeling; however requirement of equipment, scaling problem, measurement difficulties 
and operating costs are the limitations of experimental modeling. Analytical approach is 
having the advantages that it has general information in formula form; however; restriction 
to simple geometry and physics and usually restricted to linear problems are the 
limitations of analytical approach. The numerical modeling has no restriction to linearity 
and complicated physics can be treated. Time evolution of flow can be obtained in 
numerical modeling. Truncation errors, boundary conditions problems and the computer 
cost are the limitation of numerical modeling. Although, the experimental wok is more 
realistic, the numerical modeling is used in this study to simulate flow in two channels 
with link. 

In this study, flow in two channels with a link is analyzed for both rigid bed and mobile 
bed channels. For the purpose, a numerical model is developed to simulate unsteady flow 



with shocks. The developed model is then applied to simulate various flood problems; 
single reach channel (rigid bed), two channels with link (rigid bed), single reach channel 
(mobile bed) and two channels with link (mobile bed). 

1.1 REVIEW OF LITERATURE 

Flood routing is described as the procedure whereby the characteristics of a flood wave at 
one location along a channel are determined from known (or assumed) data at an 
upstream location (Mahmood and Yevjevich 1975). More generally, routing refers to the 
mathematical description of water movement (in terms of flow rate or depth), from one 
location to the other. A flood is most often the result of a volume of runoff entering a 
channel with in restricted period of time, resulting in a translation wave of long wave 
length and small amplitude relative to its base length which moves downstream, changing 
in characteristics as it progresses. 

The object of flood flow computation is to determine one or more of the following 
characteristics of a flood wave: 

1. Height -the maximum surface elevation reached; the rate of rise and fall 
(subsidence) of surface elevation; and the timing involved. Such information is 
needed in the planning or design of structures across or along streams or rivers. 

2. Peak discharge -the maximum volume rate of flow is needed in the design of 
spillways, bridges, culverts, and channel sections. 

3. Total volume of flow -the volume of water involved in a flood is significant in the 
design and operation of storage facilities for flood control, irrigation, water 
supply. 
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4. The area inundated —related to maximum stage, the flood area along any stretch of 
river is of obvious interest with regard to occupation of flood plains and structure 
built on them. 

5. The timing of flood - question of interest regarding time have to do with the 
duration of flood conditions and the warning time between the initiation of flood 
conditions in the channel and the attainment of dangerous or damaging stages at a 
particular location. 

Generally all flood routing techniques require certain data on flow conditions and channel 
characteristics: 

A. Channel characteristics 

Topography of reservoir or river and neighboring valley, including contours below 
normal water level.; Storage volume in a reservoir or stream valley as a function of water 
surface elevation.; Channel bottom slope. So-', Channel cross section at boundaries of the 
reaches to be used in the model including the flood plain.; Resistance characteristics of 
the channel in terms of Manning’s n, Chezy’s C, or conveyance including some 
evaluation of the variation of resistance with stage. 

B. Flow characteristics 

Past flow records to establish stage-storage, stage-discharge, and discharge-storage 
relationship; Measurements or estimates of local overland inflows into the total length of 
channel, precipitations records, soil moisture indices, runoff characteristics, drainage 
basic areas, and other hydrologic data involving the rainfall -runoff relationship. 
Tributary flows may be measured, estimated, or (rarely) determined from a separate 
routing procedure; Measurements or estimates of the volume inflow to the upper end of 
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the total channel length as a function of time. Depending upon the technique used, stage 
verses time data, discharge verses time, or a discharge verses stage curve may be required 
as a boundary condition or at other location.; Past record of stage-time and/or discharge 
time data. 

1.1.1 Channel Networks 

An alignment of channels one the ground is channel network. In open channel flows, 
there are different types of channel networks as shown in figure 1.1. Choi and Molinas 
(1993) presented a simulation algorithm for one dimensional unsteady flow routing 
through channel networks. Khan (1993) presented Muskingum flood routing model for 
multiple tributaries, in which the Muskingum storage equation was redefined to include 
the tributary flows as independent inflows to the channel system by increasing the 
number of x parameters and the routing equation is accordingly modified to contain 
tributary flows as individual inflows through an increased number of coefficients. Nguyen 
and Kawano (1995) gave an alternative double-sweeping method applicable to a range of 
network types. Ji (1998) introduced a different type of discretization scheme for the 
complete unsteady flow equations and an algorithm to solve these equations for 
sewer/channel networks with both free surface and surcharged conditions. Ping and 
Xiaofang (1999) gave a hydraulic flood routing method for multi-Branch Rivers by 
means of the imaginary channel length using the double-sweeping method. Wu et al. 
(2004) proposed a one-dimensional model to simulate the nonequilibrium transport of 
nonuniform total load under unsteady flow conditions in dendritic channel networks with 
hydraulic structures. The equations of sediment transport, bed changes, and bed-material 
sorting were solved in a coupling procedure with a direct solution technique, while still 
the model was decoupled fi-om the flow model. This coupled model for sediment 
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calculation was more stable and was less likely to produce negative values for bed- 
material gradation than the traditionally fully decoupled models. 




(c) 

Figure 1.1: Definition sketch; Channel network 

It is seen that most of these schemes are implicit. In an implicit scheme unknown appears 
in the difference equation in a manner that gives evolution in terms of unknown 
quantities. It is well known that, although the implicit scheme will yield inherently stable 
algorithms for the solution of an open-channel network problem, a large number of 
algebraic equations may have to be solved simultaneously at each time step; as well, a 
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large number of iterations may be required in order to solve the system at a desired 
accuracy. Moreover, one still needs to use time steps of the order of the spatial grid size 
even in an implicit scheme to achieve a uniform accuracy in both space and time 
discretizations. This is exactly the order of the time step used in an explicit scheme (Aral 
et all 998). 

To overcome this limitation of implicit schemes, explicit schemes are used in channel 
networks. In the explicit schemes; only one unknown appears in the difference equation 
in a manner that permits evolution in terms of known quantities. 

Jin and Fread (1997) presented a characteristics-based, upwind, explicit numerical 
scheme for one-dimensional (ID) unsteady flow modeling of natural rivers and 
implemented into the U.S. National Weather Service (NWS) FLDWAV model in 
combination with the original four-point implicit scheme. This explicit scheme was 
extensively tested and was compared with the implicit scheme. The study showed that the 
new explicit scheme provided improved versatility and accuracy in some situations, such 
as particularly large dam-break waves and other unsteady flows with near critical mixed- 
flow regimes. A technique for implicit-explicit multiple routing to incorporate the 
advantages of using both the schemes within an application of the FLDWAV model. Aral 
et al. (1998) presented use of nonoscillatory second order relaxation scheme to simulate 
the flow in open channel networks. Junction flow conditions in the channel network were 
solved explicitly using the continuity principle at each junction and the characteristic 
equations. In that study, flow in two hypothetical channel networks was presented and the 
model gave satisfactory comparison with widely used Preissmann four-point implicit 


scheme. 
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In most of the implicit and explicit schemes, junction conditions are solved by continuity 
equations and the assumption of equal flow depths at the three or more cross sections of 
junction. It is shown by Garcia-Navarro (1993) that this is a valid approximation only for 
few simplified cases. He presented method for surges through an open channel junction 
and treated the surges through the junction as a hydraulic jump. Assumptions for the 
method were- the bottom slopes are all zero and the friction is negligible; the channels are 
of uniform rectangular cross section; only the junction in which two streams combine in a 
single channel has been considered; the flow is parallel to the channel walls immediately 
above and below the junction; and the widths of the channels verify: b 3 = bi+b 2 . It may be 
noted that in the present study, the assumption of equal depths are used. 

1.1.2 Governing Equations 

To show any physical problem in mathematical form, governing equations are necessary. 
Different types of governing equations may be used in open channel flows. A three- 
dimensional fluid flow is represented by the Navier-Stokes equation. Newtonian fluid is 
the only assumption used in this equation. But the solution of Navier-Stokes equation is 
not simple: Generally St. Venant equation is used as the governing equation for flow in 
open channels. Assumptions used in this equation are present in chapter II. Further 
simplification in momentum term of St. Venant equation results in the 
kinematic/diffusion equations. 

Although, Open channel flows are not always represented by one-dimensional e.g. curved 
flows, flow at channel expansion, junction flow etc. However, flow is assumed to be one- 
dimensional in the present study. The sediment transport capacity model consists of the 
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mass conservation or continuity equations for the water-sediment mixture and bed 
material (sediment) as well as the momentum conservation equation for the mixture (Cao 
et al. 2004). By taking the density of water equal to the density of sediment mixture, the 
continuity and momentum equation for water is not changed as given by St. Venant 
Equations. Bed elevation change is ignored in the continuity equation for water. Sediment 
continuity equation is more simplified by eliminating change of rate of sediment density 
in mixture (Saiedi 1997). For the present study the simplified form of sediment continuity 
equation (Saiedi 1997) is used. 

1.1.3 Numerical Methods 

The first attempt to numerically simulate the one-dimensional St. Venant equations for 
flood routing problems was made by Stoker (1953a) and Isaacson (1954). In 1961, 
Preissmann provided the famous four-point box scheme. Since than, a variety of 
numerical techniques have been used to solve flood routing problems. Many of these 
standard numerical models are discussed or described by Chaudhry (1993). After 1990, a 
variety of high-order shock-capturing schemes have been used to solve St. Venant 
equations. Many of these shock capturing models are discussed or described in papers or 
textbooks (Toro 2001; Zoppou and Roberts 2003). In the present study explicit second 
order accurate shock capturing TVD scheme in finite-difference approach is used. By 
changing parameters related to TVD scheme this scheme can be changed to essentially 
non-oscillatory (ENO) schemes. At the boundaries and junctions the method of 
characteristic is used. 

The functions of continuous argument which describe the state of flow are replaced by 
functions defined in a finite number of grid points within the considered domain. The 
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derivatives are then replaced by divided difference. Thus, the differential equations, i.e., 
the laws describing the evolution of continuum, are replaced by algebraic finite difference 
relationship. The computational grid is a finite set of points, sharing the same domain in 
the {x, t) plane. 

Most traditional models fail to predict the transcritical flows because of numerical 
instability. Therefore, High order shock capturing schemes have been extended from the 
non-linear scalar case to one-dimensional systems, based on the use of Riemann solvers. 
There is various type of higher-order shock capturing schemes. 


1.1.4 TVD Schemes 

The concept of TVD scheme was introduced by Harten (1983). The theory which is used 
to establish a total variation diminishing scheme is based on scalar equations. A measure 

of the oscillations in a solution is provided by the total variations 

- 1-00 

TV(U") = X I ^j" I values Uj at the time level n. Oscillation are avoided 

7 =-“ 

if the total variations do not increase with time. A general form of an explicit scheme can 
be written as 







The numerical flux for the general scheme can be written in the form 

■^± 1/2 •^<■± 1/2 A ± 1/2 ) 

Where R i±i /2 is the right-eigenvector matrix corresponding to the approximated Jacobian 
J, and Di±i /2 is the scheme-dependent vector function. Element of vector Di±i /2 is 
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different for different type of scheme (Delis and Skeels 1998). Also the element of vector 
Dj±i /2 have limiter function applied to the scheme. The purpose of the limiter function is 
to supply artificial dissipation when there is a discontinuity or a strong gradient, while 
adding very little or no dissipation in regions of smooth variation. Some of the limiters 
(Zoppou and Roberts 2003) are given in table 1.1 


Table 1.1: Example of limiters applicable to solved second- or higher-order numerical 


schemes for solving conservation laws 


Limiter 


Minmod 


Superbee 

max 1^(0, min(l, 2y0"), min(2, 

+ max [(0, min(l, 2;5'^), min(2, 

MUSCL 

max|o, ,{I3* + 

Van Leer 

\+\P^\ \+\p-\ 

Van Albada 

1+irf 

Monotonic 

p^+\p^\ ^p-+\p-\ ^ 

l+\P^\ 1+|;0-| 

Here the , 

^ ^./>3/2 ^^i+1/2 ^y+i/2 “ ^y+i '~^j 


> Ol 


Although application of TVD to nonlinear systems does ensure that the solution will be 
total variation diminishing, however TVD schemes used to one- dimensional open 
channel flow equations with source terms, successfully. Garcia-Navarro et al. (1992) used 
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the TVD-MacCormac scheme for simulating 1-D open channel flow. Delis and Skeels 
(1998) used the TVD scheme for open channel flows with source terms. The performance 
of second- and third- order-accurate TVD schemes was investigated for the computation 
of free surface flows, in predicting dam-break and extreme flow conditions created by the 
river bed topography. Tseng et al. (2001) presented two high-resolution, shock-capturing 
schemes for the simulation of ID, rapidly varied open-channel flows. The scheme 
presented incorporates the method of characteristics to deal with the unsteady boundary 
conditions. Also, the Strang-type splitting operator was used to include the effects of 
bottom slope and friction terms. The results of dynamic flood routing and steady routing 
was compared to demonstrate the risk of using steady routing for flood mitigation. 
Burguete and Garcia-Navarro (2001) presented the new technique in which source term 
was included in the flux limiter functions to get a complete second-order compact 
scheme. Sanders (2001) presented high-resolution and non-oscillatory solution of the St. 
Venant equations in non-rectangular and non-prismatic channels. Tseng (2003) presented 
four different approaches to discretize the source terms for the simulation of one- 
dimensional open channel flows with rapidly varied bottom topography using TVD- 
MacCormack-type predictor-corrector scheme. Compared with other high-resolution 
shock-capturing schemes, MacCormack-type predictor-corrector method is easy to 
implement and does not presents any additional difficulties in dealing with the source 
terms. Cao (2004) used TVD scheme for computational dam-break over erodible 
sediment bed. This paper presented that in early stage of the dam-break, a hydraulic jump 
is formed around the dam site due to rapid bed erosion, which attenuates progressively as 
it propagates upstream and eventually disappears. It was also presented that the free 
surface profiles and hydrographs was greatly modified by bed mobility, which has 
considerable implications for flood prediction. 
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1.1.5 ENO Schemes 

Shu and Osher (1988) and Nessyahu and Tadmor (1991) proposed the non oscillatory 
schemes. As presented by Nujic (1995) the simplified version of ENO schemes are 
presented here. In the ENO scheme the numerical flux is splitted up into two parts 

■^y+l/2 ~ ■^j+I/2 ■*■■^7+1/2 

For second-order accuracy numerical fluxes are defined as 

F;,y,=F;+o.5SF; 

•^7+1/2 ~-^/+i ~0.5SFj^^ 

Where 


F; = 0.5(F, +aUj),F- = 0.5(F, -aUj) 

Here a represents some positive coefficient, and it is required that 

a > max I /I, I i = l, ,m 

j 7=1, 

m, number of unknowns and N total number of grid points. 

sf ; = min mod(F,;, - f ;, f ;- f ;_,) 


The quantities 


SFJ^^ = min mod(F/^, - FJ , FJ^^ - FJ^, ) 


And each of the fluxes is approximated up to the required order of accuracy. 

Tseng (1999) presented the comparison between TVD and ENO schemes. TVD and 
ENO schemes were shown in a single derivation by using three different parameters. 
Tseng et al. (2001) presented ENO scheme for the simulation of ID, rapidly varied open- 
channel flows. Other high order shock capturing schemes such as upwind schemes (Ying 
et al. 2004), Relaxation scheme (Jin and Xin 1995 and Aral et al. 1998), flux vector 
splitting scheme (Jin and Fread 1997), the modified Godunov method (Savic and Holly 
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1993), the space-time conservation schemes (Chang 1995) and used by Molls and Molls 
(1998), Lagrangian discrete parcel method (Wang and Shen 1999) and Petrove-Galerking 
finite elemant scheme (Yang et al. 1993; Khan 2000) were applied to open-channel flows. 

Toro (2001) provided an indepth review of TVD schemes and approximate Reimann 
solvers for the shallow water equations. Zoppou and Roberts (2003) summarized 20 
different high resolution shock capturing schemes for dam-break flow and concluded that 
verity of shock-capturing numerical schemes were efficient, accurate, robust, and suitable 
for solving the shallow water equations when discontinuities were encountered in the 
problem. It should be noted that the TVD scheme has second order accuracy, but the 
characteristic equations which are employed at the boundaries and junctions are only first 
order. MacCormack heuristically showed that if the order of accuracy of the boundary 
conditions is an order less than that of interior points, than the overall accuracy of the 
computed results is not impaired, although other researchers questions the validity of this 
statement (Fennema and Chaudhry 1986). 

1.1.6 Sediment Modeling 

In case of water flow in a channel with mobile bed it is require to solve the sediment 
continuity equation with the water continuity and water momentum conservation 
equation. There are three different approaches for solving these equations and find the 
value of flow discharge, area and bed elevation to next time level. In a coupled approach, 
the physical coupling of water and sediment phases of alluvial flows is recognized and all 
equations including those related to both water and sediment are solved simultaneously. 
In the uncoupled approach strong interaction between the bed elevation change and water 
flow is eliminated. In the uncoupled approach water continuity and water momentum 
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conservation equation solved simultaneously and the sediment continuity equation is 
solved separately. If the time step for water flow equation and sediment continuity 
equation is same in the uncoupled approach, then it is known as semi-coupled model. The 
main draw back in most uncoupled and semi-coupled models are the assumption of a 
fixed bed channel when solving flow equations and lack of inclusion of sediment 
parameters in determining flow resistance. 

In many practical cases the use of uncoupled approaches may be justified, mainly because 
of different time scale in motion of sediment and water, plus the fact that our information 
on rules of sediment transport and hydraulic resistance in mobile rivers involved 
empirical relations with various degrees of approximation that largely influence the 
accuracy of the results from any uncoupled and coupled model. In modeling movable bed 
channels subjected to peaky floods or strong tidal currents, significant and rapid changes 
in bed levels, a decoupled model may lead to unacceptable errors in long term prediction 
of river changes or even to computational difficulties. Saiedi (1994) showed that if the 
relations used for sediment transport and hydraulic friction in a coupled model are the 
same as those used in another uncoupled model, the coupled model gives more realistic 
results. Gill (1983), Nastha et al. (1987), Zhang and Kahawita (1987) and Mohapatra 
(1991) used decoupled modeling. Mohapatra and Bhallamudi (1994) presented the 
analytical and numerical model for the analysis of steady-state bed level variation in 
gradually expanded open channels with movable bed. A second-order accurate, explicit 
finite-difference scheme (MacCormac scheme) with uncoupled approach had been used 
by them. Bhallamudi (1991) presented a coupled model for the flow simulation in mobile 
bed channels. In that study, the MacCormack explicit finite-difference scheme was used 
which, was second-order accurate, handle shocks and discontinuities in the solution 
without any special treatment, and allows simultaneous solution of the water and 
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between results of uncoupled a coupled model. He presented different scheme for 
discretize the sediment continuity equation and compare with coupled model. Capart and 
Young (1998) presented the formation of jump by the dam-break wave over a granular 
bed. In this study a numerical model was present with experimental results and show that 
shock formation involves a particularly strong coupling between flow free-surface 
evolution and bed morphodynamics. Cao et al. (2002) presents coupled and decoupled 
numerical modeling of flow and morphological evolution in alluvial channels. In this 
study four approaches for water and sediment flow by different simplifications both in 
water and sediment equations and concluded that it is require fully coupled model is 
require to simulate extreme flow in channels. Singh et al. (2004) presents a fully coupled 
one-dimensional alluvial river model for simulation of sharp hydraulic and bed transients 
in alluvial river. In this study the governing systems of partial differential equations are 
discretized by making use of generalized Pressmann finite difference scheme. This model 
is further extended to simulate the process of grain sorting, wash load transport and non- 
equilibrium sediment transport. Although it is shown that a coupled model is more stable, 
the semi-coupled approach is used in the present study. 

In this chapter, importance of the problem and relevant earlier works are presented. 
Governing equations are presented in chapter II and their numerical solution is presented 
in chapter III. Results of the present study are presented in chapter IV. Important 
conclusions are given in the last chapter. A brief note about interlinking project in India is 
given in appendix. 
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CHAPTER II 

GOVERNING EQUATIONS 


In the present work, flow in two channels with a link is analyzed by using a numerical 
model. With reference to the present study, governing equations are comprised of 
continuity and momentum equations for water flow and continuity equation for sediment 
flow. These equations are described in the following sections. 

2.1 WATER FLOW 

The fundamental equations of fluid dynamics are based on the following universal laws of 
conservation: Conservation of mass; Conservation of momentum; and Conservation of 
energy. However, in this study, only the conservation of mass (continuity equation) and 
the conservation of momentum laws are used to describe water flow. 


For the one-dimensional unsteady gradually varied flow in open channels, continuity and 
momentum equations (St. Venant 1871; Delis and Skeels 1998; Tseng 1999; Burguete and 
Garcia-Navarro 2001) are; 

Continuity equation'. 


dt 


2.1 


Momentum equation in longitudinal direction for water 


dQ ^ d{Q^IA + gf) 
dt dx 


gA{-^-Sf) + gh 


2.2 


In the above equations (2.1 and 2.2), t = time; a: = distance along the channel; A = wetted 
cross-sectional area; Q = volume flow rate; g = gravitational acceleration; Z = bed- 
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elevation; Sf— friction slope; and Ij and I 2 represent the hydrostatic pressure force and the 
pressure force due to width variations, respectively. 


The friction slope is calculated using the Manning equation 
. 2/^2 




-> "I 


2.3 


In which, n= Manning roughness coefficient and R, hydraulic radius, is the ratio of the 
cross sectional area. A, and the wetted perimeter, P. The hydraulic radius is equal to the 
flow depth for wide rectangular channels. Although, n is a complicated function of flow 
depth, bottom roughness, bed slope, discharge and bed forms, a value independent of these 
factors are assumed throughout the study. 

The hydrostatic pressure force (/;) is calculated using the following equation 

h{x,t) 

7, = J (h-r\yy(x,T\)dr\ 2.4 

0 

and the pressure force due to width variations (h) is calculated using 

72= J (h-v[) — ^-^dr\ 2.5 

Where, b(x, t]) = channel width at depth q; and h = total water depth. For trapezoidal cross 
sections, 7/ and h are given by 




2.6 


h 


,2,1 dB hdS,.^ 
r( + 

2 dx 3 dx 


2.7 


Where B = width of channel bottom, and Sl = side slope (vertical : horizontal) of the 


channel. 
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Where B - width of channel bottom, and Sl = side slope (vertical : horizontal) of the 
channel. 


r- 1 




Cross section at 1-1 

Figure 2.1: Definition sketch: unsteady flow in open channels 

The assumptions used in equations 2. 1-2.2 are: 

(i) Flow is one-dimensional in all channels. It may be noted that flow at the 
junction is not one-dimensional due to stream line curvature. However it is 
assumed to be one-dimensional in this study. 

(ii) Slope 9 of the channel bottom is small so that Sin6 = Tan0 = 0 and Cos0 = 1 . 

(iii) The pressure distribution at a section is hydrostatic. This is true if the vertical 
acceleration is small, i.e., if the variation of water surface with alignment is 


gradual. 
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(iv) The transient-state friction losses may be computed by using formulas for the 
steady state friction losses i.e., the roughness coefficient formulas can be used 
to find the slope of energy line. 

(v) The velocity distribution at a channel section is uniform i.e. the velocity is 
uniform over the cross section and the water level across the cross section is 
horizontal. 

(vi) The channel is straight and prismatic. 


Equations 2.1 and 2.2 can be represented in a vector form (Tseng and Chu, 2001) as 


dU 

dt 


dx 


2.8 


Where, U = 




;G = 


0 

gA(S,-S,) + gI, 


2.2 JUNCTION FLOW 

Junctions in a channel network are treated as internal boundary conditions. At the 
junction, there may be combined or separate flow. 



Figure 2.2: Junction conditions 
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In addition, energy equations between sections 1-2 and 1-3 are also to be satisfied. 






2gA^ 


2gAr 




2gA3^ 


Where h/is head loss due to friction. 

2.3 SEDIMENT FLOW 


2.10 


2.11 


Like water flow, the sediment continuity and momentum equations are required to 
represent the sediment flow in channels. Sediment continuity equation can be represented 
(Saiedi 1997) as: 


dZ 1 dq^ _ 
dt (1 — p) dx 


2.12 


Where, q,, = total sediment discharge per unit width, and p = the porosity of bed material. 


The sediment flow in a channel is a very complicated phenomenon. The physics of 
sediment flow is not well understood. Therefore, in place of momentum equation for 
sediment flow, a simple relationship between the water and sediment discharge is used in 
addition to the sediment continuity equation (Eq. 2.12), to completely represent the 
sediment flow. The flow discharge and depth, particle sizes of bed material and bed-forms 
are the factors on which the sediment discharge depends. Sediment discharge is divided 
into three parts: bed load, suspended load and wash load. The bed load is the part of 
sediment discharge which is near the bed and suspended load is the part of sediment 
discharge in form of suspended material in water. The total load is the summation of the 
bed load and suspended load. The wash load is generally due to erosion of side walls of 
channels. The bed load is the major part of the sediment discharge in a channel. Several 
equations are available for estimating bed load, suspended load, wash load and total load 
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channels. The bed load is the major part of the sediment discharge in a channel. Several 
equations are available for estimating bed load, suspended load, wash load and total load 
in tenns of water flow (Graf 1972, Garde and Rangaraju 1985). As suggested by Vanoni 
(1975) and Jansen et al. (1979) many of these equations can be represented in the 
following functional form. 

2.13 

Where, is a transport parameter and 7 is a flow parameter. The transport parameter 
depends on the sediment discharge and the grain properties, while the flow parameters 
depend on the flow properties and the grain size. 

For simplicity, the following sediment transport formula is used in this study (Gill 1983, 
Zhang and Kahawita 1987, Park and Jain 1987, Mohapatra 1991). 

^.,=c(% 2.14 

A 

Where, c and d = empirical constants that depend on the sediment properties and the bed 
forms. This equation is valid only for uniform sediment size. It may be noted that no 
distinction is made between the bed load and the suspended load and equation 2.14 is 
assumed to represent the total sediment discharge. The aim of this study is to analyze the 
flow in two channels with a link. Therefore, the use of a simplified relationship is 
justified. For real field problems, an appropriate sediment transport equation in place of 
equation 2.14 should be used. 
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CHAPTER III 

NUMERICAL SCHEME 

The continuity and momentum equation for water, and the continuity equation for 
sediment flow are the governing equations used in this study. The governing equations 
(2.1 and 2.2) constitute a set of non-linear hyperbolic partial differential equations. In 
these equations, x and t are the independent parameters; discharge (Q), cross sectional 
area {A), and bed-elevation (Z) are the dependent parameters. The dependent parameters 
are known for all x at the initial time level and these are to be calculated for the next time 
level. Analytical solutions for these equations are available only for idealized cases. 
Therefore, they are solved numerically. In this study, a finite-difference scheme is used 
for solve these equations numerically. The functions of continuous argument which 
describe the state of flow are replaced by functions defined in a finite number of grid 
points within the considered domain. The derivatives are then replaced by divided 
difference. Thus, the differential equations, i.e., the laws describing the evolution of 
continuum, are replaced by algebraic finite difference relationship. The computational 
grid is a finite set of points, sharing the same domain in the {x, t) plane. A typical 
computational grid for simple one-dimensional problem is shown in Fig. 3.1. In this 
study, the governing equations are solved in a semi-coupled approach, i.e. the continuity 
and momentum equation for water flow are solved simultaneously and discharge and 
cross-sectional area are used to solve sediment continuity equation. 
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Figure 3.1: Finite-difference grid 


In the present study, second order accurate Total Variation Diminishing (TVD) scheme is 
used to solve the governing equations 2.1 and 2.2. TVD is a high order shock capturing 
scheme which, numerically captures discontinuities with sharp comers and avoid 
unrealistic oscillations. For the governing equation 2.8, this algorithm ensures that the 
total variation does not increase with time, that is: 

/ / 

Where, , - U" . The solution can be second- or third-order accurate in non 

critical sections, but switches to first order at extreme. 


To explain the numerical method, any dependent parameter is expressed asF" , where i 

and n represent the indices for node number and time level. For example, Q° is the initial 

value of discharge at node number 3. As mentioned above, the dependent parameters, at 
initial time level, are known and it is required to obtain their values at subsequent time 
levels. The solution strategy is shown in Fig 3.2 and details are provided in the following 


sections. 
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Figure 3.2: Flow chart for numerical model 
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3.1 INPUT PARAMETERS 

The input parameters for the model are cross-sectional geometry. Manning roughness 
coefficient (n), length of channel L, number of divisions K, acceleration due to gravity 
(sX porosity of bed material, empirical constants a, c and d (equation 2.14), Courant 
number (Cn) and time level at which results are required. In addition, the boundary 
conditions at the inflow and outflow section are given as input. 

3.2 GRID PARAMETERS 

Ax is given by 

Ax = ^ 3.1 

K 

Referring to Fig. 3.1, i,T,in= 1 and in,ax=K+l. 

3.3 INITIAL CONDITION 

Initially (t=0) the dependent parameters (Q, A, Z) are either known or are derived from 
the input parameters. 


3.4 STABILITY 

For the stability of the present numerical scheme, Courant-Friedrich-Lewy condition is 
satisfied: 


At, = C * 

^ i n 


Ax 




3.2 


Where celerity, 



. Time step. At, is 


chosen as the minimum of all Ati. 
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3.5 UNSTEADY FLOW CALCULATION 

3.5.1 Water Continuity and Water Momentum Equation 

Total Variation Diminishing (TVD) scheme is used to solve the governing equation 2.8. 
The homogeneous part of the system is hyperbolic and is responsible for most of the 
difficulties when it is numerically integrated. The non-linearity of the flux vector F (Eq. 
2.8) may lead to spontaneous oscillations, which may have real physical meaning. 
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With non-zero source term, Strang-type splitting operator (Strang 1968) is employed to 
Eq. 2.8 (Seng et al. 2001). 


U"*' = L 

I .V 

Where, 




V ^ y 


L{At)L^. 




J 


U" 


3.7 


L,{At)U” =U” +AtG^ + 
and 


(Af)^ f dOX 
2 



3.8 


The numerical flux vector > ^nd the source term averages, G" depends on various 
schemes. The numerical flux for the general scheme can be written in the form 

■*■■^±1 ■*'-^;±1/2A±1/2) 3.10 

Where R i±i /2 is the right-eigenvector matrix corresponding to the approximated Jacobian 
J , and Di±i /2 is the scheme-dependent vector function. 


For TVD2 scheme, the component of Di+ 1/2 is given as 

A+r/2 ~ ^(^/+I/ 2 )(^|+I )~ 9 (\+ 1/2 ■*■71+1/2) ®i+l/ 2 ’ J = 3.11 

Where, 

i«,dAC/] 3.12 

<j(z)-iL(z)- Azd 3.13 

2L Ax _ 

~ ^[®i+l/2’^i-l/2] 



1 _ j'^(^,+i/2)(^M if a!^ip5iO] 

U+Ml 1 

0 otherwise 


The limiter function is given as 

mia, Z)) = ^ I if sgn a = sgn b = s] 

[ 0 otherwise j 

The entropy correction function is given as 


(p(z) = - 


e <z 


(z^+s^)/ 28 if|z|<s 


Here, e is a small positive number. 
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3.15 


3.16 


3.17 


3.5.2 Boundary Conditions for Q and A 

The model presented in Eq. 3.7 gives the values of and at node 3 to K-1 in a 
channel. Here, K+1 is the last node of channel reach. For the channel networks, at the 
junctions Q and A can not be found by TVD scheme; therefore boundary conditions are 
required. 

Open boundaries 

The node number 1-2 and K-\ and K are open boundaries for a channel. For node number 
1, flow hydrograph or stage hydrograph is the physical boundary condition. For down 
stream end flow hydrograph, stage hydrograph, normal depth, or rating curve (discharge 
Vs stage) can be given as physical boundary condition. Other than the physical boundary 
conditions, characteristic method (Chaudhary 1993, Garcia-Navarro and Saviron 1992, 
Aral et al. 1 998) or some approximation is used to find Q and A at the boundaries (node 
nos. 1 and K+1). 
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In method of characteristic (MOC), positive and negative characteristic equations are 
written as; 


es" 

and 

a"" =C„+C„, 

In which, 

C'r=a".+C,-,y",+«W,-irV,At 

^ H = Qm ~^o/+i>’,+i At 


3.18 

3.19 

3.20 

3.21 

3.22 


For node numbers 2 and K, two characteristics equations (one positive and one negative) 
are used to find the values of Q""^' and A"""’. Alternatively, for these nodes condition of 
node no. 1 and K+l may be used. 


Junction 

By referring Fig. 2.2 and applying equation 2.10 

q;’"'=q^*'+q f 

Assumption of equal depths at 1, 2 and 3 gives: 

h;'"'=hf' =hr' 

Eqs. 3.18 and 3.19 gives 
Q;’"'=f(h;'"') 

Qf =f(hr') 


3.23 


3.24 


3.25 (a) 
3.25 (b) 


3.25 (c) 
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By substituting 3.25 (a), 3.25 (b) and 3.25 (c) in Eq. 3.23 gives a common value h and 

then by equations 3.25 (a), 3.25 (b) and 3.25 (c) gives three values of Q vales at the 

junction. 


3.5.3 Sediment Continuity Equation 

The third dependent variable, Z, is determined at the next time level by solving equation 
2.12. As mentioned above, semi-coupled model is used to find the bed level changes, i.e. 
sediment discharge is found by discharge and depth values and is used to estimate the bed 
level change in the channel. The following finite difference approximations are used to 
solve the sediment continuity equation (2.12). 


”1 n'fiM yn 

dt _i Ai 


3.26 


dx 






n+l 


2Ax 


3.27 


Substitution of equations 3.26 and 3.27 in 
at the unknown time level, 


f -cv'i 

f 

1 2*i\-p) J 

vAxJ 


equation 2.12, gives an explicit equation for Z 
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3.5.4 Boundary Condition for Bed Elevation 

Equation 3.28 is applicable for points 2 to K. Therefore, Z at node number 1 and K+1 is 
determined by boundary conditions. In addition for the junctions, it is required to modify 
the equation. The bed boundary condition is supplied in the form of the known sediment 
inflow rate as a function of time. In order to express sediment inflow in terms of Z values, 
the following equation is used (Singh et al. 2004): 
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3.29 


Where, q^.-^ is the sediment inflow at the unknown time level. Bed level at the 


downstream end is obtained by extrapolation from the interior points as; 

Z Aj+1 /-> 'T'W-bl o-rt+l ^ ^ ^ 

jf+i=2Zj^. -Z^_, 3.30 

After determining Q"^', A"^', Z"'*’' for all nodes, time level is increased by At. This new 
time level is compared with the required time level. If the required time level is not 
reached then reassign the current values of the dependent parameters as the known time 
level and repeat the earlier steps (Fig. 3.2); else, the computation is stopped depending on 
the requirement, appropriate printing of the dependent variables is performed. 




CHAPTER IV 

Results and Discussion 
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In the present study, a one-dimensional model using the total variation diminishing 
scheme (TVD) has been applied to study the flood propagation for two channels with a 
link. Three different cases, as given below, are studied to meet the objectives 

(i) Single channel (rigid bed), 

(ii) Two channels with a link (rigid bed), 

(iii) Single channel (mobile bed), 

Results for these cases are presented below. 


4.1 SINGLE CHANNEL (RIGID BED) 

Before applying the present model to study flood routing in a single channel, the 
developed model is validated against the analytical results for dam-break flows (DBF). 
Because, DBF are with shocks, which provides a check for the shock capturing 
capabilities of the present methodology (TVD scheme). 

4.1.1 Model Validation 

The developed mathematical model (as discussed in chapter III) is used to simulate 
Stoker’s solution for dam-break flow (DBF). Stoker (1957) presented analytical results of 
dam-break flow by using following assumptions: (i) Channel bed is horizontal; (ii) Both 
reservoir and the channel are infinitely long; (iii) Channel bed and walls are frictionless; 

(iv) the dam failure is instantaneous; (v) prior to the dam failure, water in both reservoir 
and the channel is static (Fig. 4.1). 
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(a) Before the dam-break 



(b) After the dam-break 

Figure 4.1: Definition sketch: Dam-Break flow 

Different input parameters used for the purpose are: X=1000m; So =0.0; «=0.0; hu=10.0m; 
hd=5.0m and time after dam failure=20s. Numerical parameters such as step size. Ax and 
Courant number, Cn are varied. It may be noted that these set of input parameters refer to 
the assumptions used in Stoker’s solution. The results indicate four well defined zones in 
the surface profile after the dam-break. Effects of the various parameters are discussed in 
the following subsections. 
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4.1. 1.1 Effect of grid size (AX) 

It IS traditional to verify any numerical results for grid independence. The effect of grid 
size, on the computation of surface profile is presented in Fig. 4.2. Three different values 
ot the step size (Ax=5m, 2m, Im) is used to find the water surface profile. In this figure, 
the surface profile is indicated only at the vicinity of the shock. 



675 685 695 705 

DISTANCE (m) 

Figure 4.2: Effect of grid size (Cn= 1 .0) 

As shown in figure 4.2, the best match with the analytical results is predicted by Ax = Im. 
Other values of Ax result in smearing of the sharp wave front. Further decrease in grid 
size does not improve the results significantly. 

4.1.1.2 Effect ofC„ 

Fig 4.3 shows the effect of Cn on the computation of the water surface profile due to 
DBF. 
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Figure 4,3: Effect of Cn (Ax = Im) 

As shown in Fig. 4.3, a higher value of Cn improves the capturing of the wave front. It 
may be seen from this figure that even for a value of Cn=1.0 the present model is 
executable. Therefore, the present model is more efficient compared to other explicit 
schemes. 

4. 1.1. 3 Effect of depth ratio 

Depth ratio, r is defined as the ratio of downstream to upstream depth {r=hd/hu). In figure 
4.4, computed surface profiles by using the present numerical model for three different 
depth ratios (r=0.5, 0.2 and 0.1) are presented. In all the cases, four zones of Stoker’s 
solution are simulated satisfactorily. As expected the present numerical model captures 
the sharp wave front very accurately. It is important to note that the present model cannot 
simulate the DBF for r=0, i.e. a dry bed downstream of dam. 





Figure 4.4: Surface profile of DBF due to different depth ratios (Ax - Im, Cn 1 .0, 

time=20s) 
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4.1.2 Model Application 

1 he numerical model presented in the previous chapter and validated in the previous 
section is used to simulate the flood propagation in a single reach channel. For the single 
reach channel (rigid bed), results are compared with those obtained from HEC-RAS 
model (USAGE 2005). It may be noted that HEC-RAS is a widely used computer 
package for flood analysis in open channels. 

4. 1.2.1 Steady flow in a single reach channel 

1 he present model is used to analyze a steady flow in a single reach channel. The inflow 
discharge is increased from a constant value of 22.58mVs to a higher value of 58.5m^/s in 
a very short time period (Fig. 4.5 (a)). The input parameters are L=2000m; B=10m; So= 
0.0001, Cn=1.0 and Ax=40m. A variable n is estimated by the functional relationship 
rt=0.04-(0.02*x)/2000. The initial discharge and flow depth are taken as 22.58m^/s and 
3.0m, respectively. A normal flow depth is imposed at the down stream boundary. Fig. 
4.5 (b), presents the discharge at the out-flow section; Fig. 4.5 (c)-(d) presents variation of 
the flow depth at inflow and out-flow sections, respectively. In addition, results obtained 
by using HEC-RAS are also presented in these figures. It may be noted that results are 
grid independent. As shown in Fig. 4.5, the model obtains a steady state flow condition 
for discharge 22.58mVs and then for 58.5mVs. The comparison between present model 
and HEC-RAS i.s satisfactory. 
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4.1.2.2 Unsteady flow in a single reach channel 

The present model is used to analyze a flood routing in single reach channel in this 
subsection. Different input data used are, L=1 0000m; B=100m; n=0.02; So=0.0001. 
Inflow discharge is 300mVs up to 3600s, increases to 900mVs in 3600s and then 
decreases to initial discharge in 7200s and remain constant (Fig. 4.6 (a)). For simplicity, 
inflow hydrograph with a triangular shape is used. A normal flow depth is imposed at the 
down stream boundary. Fig. 4.6 (b) presents depth variation with time at the inflow 
section (obtained by MOC); Fig. 4.6 (c) presents discharge variation at out-flow section 
and 4.6 (d) presents depth variation with time at the out-flow section. In addition, results 
obtained by using HEC-RAS are also presented in these figures. It may be noted that 
results are grid independent. In Fig. 4.6 (b), there is a difference between depths near the 
peak value at inflow section. In Fig. 4.6 (c) and (d) there is a fluctuation at the start of 
rising limb for HEC-RAS model. However, the present model does not have such 
fluctuation. Probably, it is due to a better shock capturing by the present model. From Fig. 
4.6 (a) and Fig. 4.6 (c), there is a decrease in peak discharge fi-om 900m7s to 625.10m^/s 
and it is also to be noticed that time to peak at out-flow section is 10807s in comparison 
of 7200s at inflow section. The overall comparison between the present model and HEC- 
RAS is satisfactory. Table 4. 1 and 4.2 shows the effect of n and effect of bed slope on the 
flood peak and time to peak at the out-flow section. As n is increased, peak is reduced and 
as bed slope is increased the peak is increased. 



WATER DEPTH (ml DISCHARGE (cumecs) 
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(a) Inflow Hydrograph 



(b) Depths at the inflow section 


Figure 4.6: Results for unsteady flow in single reach channel 
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Table 4.1 Effect of non flood peak (slope=0.0001) 


Manning’s roughness 

coefficient 

Maximum discharge at 

out-flow section 

(m^/s) , 

Time to peak (s) 

0.02 

625.1 

10807 

0.03 

580.27 

11473 

0.04 , 

548.43 

11896 

0.05 

524.73 

12219 


Tabic 4.2 Effect of bed slope on flood peak (n =0.02) 


Bed slope 

Maximum discharge at 

out-flow section 

(m^/s) 

Time to peak (s) 

0.001 

827.95 

9332 

0.0005 

769.94 

9874 

0.0001 

625.1 

10807 


4.2 FLOW IN TWO CHANNELS WITH LINK (RIGID BED) 

Two channels with a link are analyzed in this section (Fig. 4.7). This type of geometry is 
not common in natural channels. However, it may occur in the following engineering 
applications as; River networking; Irrigation canal networking; Flow diversion; and 
Channel splitting. 
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Figure 4.7 Definition sketch: two channels with link 

4.2,1 Model Validation 

In the absence of any analytical results, the present model is validated against results 
obtained by using HEC-RAS. Referring to figure 4.7, the input parameters for a 
hypothetical layout are: Ll=L3=L5=4000m; L2=L4=6000m; Bl=B2=100m; B3=20m; 
B4=B5=80m; nl=n2=n3=n4=n5=0.02; Sol=So2=0.0005; So4=So5=0.0005; Zl=100.0m; 
Z4=99.0m; Ax=l 00.0m and Cn=0.8. The constant inflow discharge in channel I and IV is 
lOOOmVs and SOOm’Vs respectively. The initial depth value used is 3.5m at every location. 
The results as computed from the present model are presented in table 4.3 in which the 
nomenclature is used as that discussed in chapter III. In addition, results obtained by 
using HEC-RAS are also presented in the same table and the comparison is satisfactory. 
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Table 4.3 Summary of results for steady flow in two channels with a link 


Flow parameters in different channels 

Present Model 

HEC-RAS 

h\1 

3.79 m 

3.88 m 

hr^^,=h21=h2>\ 

3.52 m 

3.51 m 

hT 

3.56 m 

3.51 m 

/z4" 

3.76 m 

3.71 m 

M'i., =«■„ 

3.29 m 

3.30 m 


4.34 m 

4.30 m 

!23"=!23L 

133.08 mVs 

134.12 m^/s 

01'; 

1000.00 m^ /s 

lOOO.OOm^/s 

04; 

500.00 m^ /s 

500.00 m" /s 

02;>, 

866.98 m" /s 

867.80 m" /s 

05;>, 

633.60 m" /s 

634.12 m"/s 


As shown in table 4.3, the flow depths at various points in the channels and link are 
simulated satisfactorily. A detailed study showing the effects of link parameters on the 
flow diversion in channel 3 is presented in Table 4.4. Referring to Fig. 4.7 
Ll=L5=2000m; L2=L4=5000m; Bl=B2=100m; B4=B5=60m; nl=n2=0.03; n4=n5=0.02; 
Sol= S()2=0.0002; So3=0.0002; So4= So5=0.0001; 0i=lOOOm^/s; g4,=500 mVs; 
Ax= 100.0m and C„=0.8. As shown in table 4.4 the discharge in link is increased with 
increase length and decrease with increase in n. 
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Table 4.4 Effect of link properties on the flow in link 


L3 (m) 

W3(m) 

n3 

Discharge in link (mVs) 

1000 

20 

0.02 

81.56 

1000 

20 

0.025 

80.31 

1000 

20 

0.03 

78.87 

1000 

20 

0.035 

77.3 

1000 

20 

0.04 

75.62 

1500 

20 

0.03 

84.14 

2000 

20 

0.03 

88.83 

2500 

20 

0.03 

92.88 

3000 

20 

0.03 

96.67 

4000 

20 

0.03 

102.86 

5000 

20 

0.03 

112.36 


4.2.2 Unsteady Flow in Two Channels with Link (Rigid Bed) 

Eqs. 3.23 and 3.24 are used to solve the flow distribution at a junction. The direction and 
the quantity of flow in the link are depending on the various channel parameters and flow 
parameters. In addition, the effect of link on the flood peak depends on the channel 
parameters and flow parameters. The same geometrical parameters, as in subsection 4.2.1 
are used for these cases. Discharge and flow depth at the downstream section are obtained 
from the characteristic equations and the normal flow condition at this section in both 
channel II and channel V. 


Case 1 

For case 1, inflow hydrographs for channel I and IV are given in Fig. 4.8 (a). In fig 4.8 
(b), out-flow from channel II is given. In figure 4.8 (c), out-flow from channel V is given. 
In Fig. 4.8 (b) and 4.8 (c) initial parts show discharge/depth variation due to arbitrary 
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Figure 4.8: Unsteady flow in 2 channels with a link (rigid bed) (Casel) 

(contd.) 


Case 2 

The inflow hydrographs are shown in Fig. 4.9 (a). Fig. 4.9 (b) and Fig 4.9 (c) shows the 
outflows in the channel II and V respectively. As shown in Fig. 4.9 (b) outflow discharge 
for channel II is increased from 855.0m^/s to 1189.70mVs. Fig. 4.9 (c) shows that the 
peak of outflow in channel V decreased for channels with link. Fig. 4.9 (d) shows the 
flow in link. In link, flow from channel I to channel V is taken as positive. 
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Figure 4.9: Unsteady Flow in Two Channels with Link (Rigid Bed) 
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Figure 4.9: Unsteady Flow in Two Channels with Link (Rigid Bed) (Case 2) (contd.) 
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4.3 FLOW IN SINGLE REACH (MOBILE BED) 

Most ot the natural channels are with mobile beds. Many of the man made structures are 
not properly working or achieving the goal due to lack of information about the sediment 
transport in channels. Sediment discharge is present due to movable bed; the erosion of 
side walls or a source of sediment load. As discussed in chapter II, only the bed load is 
taken in to account in this study. Due to the mobile bed, the flow conditions in channel 
are changed. As discussed in the previous chapter, the semi-coupled model is used to 
analyze the flow in a reach with mobile bed. The effect of movable bed on the flow 
conditions are presented in the following sections. 



Figure 4.10: Definition sketch: single reach channel (mobile bed) 

4.3.1 Model Validation 

There is equilibrium for sediment discharge and water flow in a channel. At the inflow 
section, if sediment inflow is equal to sediment discharge capacity there is no change in 
bed elevation at inflow section. The sediment over loading is the condition for which the 
sediment input is more then the sediment discharge capacity (Equation 2. 14). The bed 
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elevation is changed for the sediment overloading. The performance of the present 
numerical results for sediment overloading is validated against previous experimental 
(Soni et al. 1980) and numerical results (Bhallamudi 1991). Referring to the figure 4.10, 
the input parameters are: L=50m; bed width B=0.2; n=0.02; constant discharge 
inflow=0.004m /s; initial bed slope=0.00356; porosity j 3 = 0.4 and initial depth=0.05m. 
Referring to equation 2.14, empirical constants c = 0.00145 and d=5.0. The grid size, 
Ax=1.0m and Cn =0.8. These parameters refer to the input values used by Bhallamudi 
(1991). Bed elevation at the upstream point is estimated by equation 3.27. Discharge and 
How depth at the downstream section are obtained from the characteristic equation and 
the normal flow condition at this section. The bed elevation at the downstream point is 
determined by extrapolation from the computed values of bed elevation at interior points. 
Fig 4.1 1 (a) and 4.1 1 (b) show the transient bed elevation and water surface, respectively, 
after 40min of over-loading. It is observed that there is aggradation, up to a distance 
15.()m. As shown in figure 4.11 (a) and 4.11 (b) the comparison of bed elevation and 
water surface is satisfactory. For two channels with a link on mobile bed present model 
does not give satisfactory results. 
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CHAPTER V 

CONCLUSIONS AND RECOMMENDATIONS FOR FUTURE WORK 


Flood wave propagation in two channels with a link was studied in the present work. A 
numerical model was developed for the purpose by using the continuity and momentum 
equations for water flow and continuity equation for sediment flow. The water flow 
equations were solved numerically by a TVD scheme, where as the sediment flow 
equation was solved in a semi-coupled manner by a second-order accurate finite- 
diffcrence scheme. Appropriate boundary conditions for water and sediment flow were 
imposed at the inflow, outflow and junctions. 

Some of the important conclusions of the present study are listed below. 

5.1 CONCLUSIONS 

The developed numerical model was applied to study the flood wave propagation in 

A. Single reach channel with rigid bed; 

B. Single reach channel with mobile bed; 

C. Two linked channels with rigid bed. 

Single reach channel with rigid bed: 

1 . Strong shocks could be captured by the present model. 

2. The maximum discharge at the outflow section decreased with increase in channel 
bed roughness. 

3. The flood wave propagation was faster in case of higher bed slopes. 
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Single reach channel with mobile bed: 

1. Present model could simulate the bed and water surface profile for sediment 
overloading. 

2. Wave like phenomenon was observed for the bed and water surface profiles as 
compared to earlier experimental and numerical results. 

3. This might be due to the semi-coupled approach used in the numerical model. 

Two linked channels with a rigid bed: 

1 . In case ot steady flow, flow divisions depended on geometries, flow conditions and 
bed characteristics of the channels and the link. 

2. Assuming the link characteristics as variables, most important parameter 
influencing the results was length of the link. 

3. In case ol' unsteady flow, maximum outflow in a channel was decreased to 15% 
when there was no link. However, the outflow was decreased to 30% in the 
presence of a link. 

4. Reverse flow in link could be simulated by the present model. 

5.2 SCOPE FOR FUTURE WORK 

The present numerical model may be further extended for the followings: 

1 . In the present study, flow in two linked channels on rigid bed was simulated. This 
may be extended to two linked channels on mobile bed. 

2. Sediment transport equation used in the present model may be replaced by 
appropriate equation for field conditions. 

3. Use of improved surface gradient method for the TVD scheme may be attempted. 

4. Coupled modeling for water and sediment flow may be included. 



5 . Improved j unction conditions may be used in the present model. 

6. Present model may be modified for application to real life problems. 

7. Inverse modeling may be performed for optimal design of the link. 
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APPENDIX 


River interlinking in India 

In oui country, the rainfall is primarily orographic, associated with tropical depressions 
originating in the Arabian Sea and the Bay of Bengal. The summer monsoon accounts for 
more than 85 per cent of the precipitation. The uncertainty of occurrence of rainfall 
marked by prolonged dry spells and fluctuations in seasonal and annual rainfall is a 
serious problem for the country. Large parts of Haryana, Maharashtra, Andhra Pradesh, 
Rajasthan, Gujarat, Madhya Pradesh, Karnataka and Tamil Nadu are not only in deficit in 
rainlall but also subject to large variations, resulting in frequent droughts and causing 
immense hardship to the population and enormous loss to the nation. The water 
availability even for drinking purposes becomes critical, particularly in the summer 
months as the rivers dry up and the ground water recedes. Regional variations in the 
rainfall lead to situations when some parts of the country do not have enough water even 
for raising a single crop. On the other hand excess rainfall occurring in some parts of the 
country creates havoc due to floods. Irrigation water demand is also increased due to 
increase in population. It is also require a lot of amount water, which cannot be fulfilled by 
small-scale irrigation projects. There is a limit to harvest water in small-scale pond, 
domestic water harvesting unit etc. 

One of the most effective ways to increase the irrigation potential for increasing the food 
grain production, mitigate floods and droughts and reduce regional imbalance in the 
availability of water is the interlinking of rivers to transfer water from the surplus rivers to 
deficit areas. If we can build storage reservoirs on rivers and connect them to other parts of 
the country, regional imbalances could be reduced significantly and lot of benefits by way 
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ol tidditioiii-il inigdtion, domestic and industrial water supply, hydropower generation, 
navigational lacilities etc. would accrue. In India 30 links have been identified as 
technically ieasible and economically viable on the basis of pre-feasibility studies. These 
are; Mahanadi (Manibhadra - Godavari (Dowlaiswaram) link, Godavari (Inchampalli Low 
Dam) - Krishna link, Godavari (Inchampalli) - Krishna (Nagarjunasagar) link, Godavari 
(Polavaram) - Krishna (Vijayawada) link, Krishna (Almatti) - Pennar link, Krishna 
(Srisailam)- Pennar link, Krishna (Nagarjunasagar) - Pennar (Somasila) link, Pennar 
(Somasila) -Cauvery (Grand Anieut) link, Cauvery (Kattalsi)- Vaigai-Gundar link, Ken- 
Belwa link, Parbati-Kalisindh-Chambal link, Par-Tapti-Narmada link, Damanganga-Pinjal 
link, Bcdli-Varda link, Netravati-Hemavati link and Pamba-Achankovil-Vaippar link. 

Similarly, based on various water balance studies carried out for the Himalayan 
component, the link proposals identified for preparation of feasibility reports include the 
Manas-Sankosh-Tista-Ganga link, Kosi-Ghagra link, Ghagra-Yamuna link, Sarda-Yamuna 
link, Yamuna-Rajasthan link, Rajastan-Sabarmati link, Chunar-Sone Barrage link, Sone 
Dam - Southern Tributaries of Ganga link, Ganga-Damodar-Subemarekha link, 
Subernarekha-Mahanadi link, Kosi-Mechi link, Farakka-Sunderbans link, and Jogigopa- 
Tista-Farakka link. 

This national Perspective Plan would give additional benefits of 25 million hectares of 
iiTigation from surface waters, 10 million hectares by increased use of ground water, 
totaling to 35million hectares which will be over and above the ultimate irrigation 
potential of 140 million hectares and 34,000 MW of hydro-power apart from the benefits 
of flood control, mitigation of drought, navigation, domestic and industrial water supply, 
fisheries, recreational facilities, salinity and pollution control etc 
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